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ABSTRACT 

Aims. We estimate stellar masses of galaxies in the high redshift universe with the intention of determining the influence of newly 
available Spitzer/TRAC infrared data on the analysis. Based on the results, we probe the mass assembly history of the universe. 
Methods. We use the GOODS-MUSIC catalog, which provides multiband photometry from the U-filter to the 8 fim Spitzer band for 
almost 15 000 galaxies with either spectroscopic (for a; 7 % of the sample) or photometric redshifts, and apply a standard model fitting 
technique to estimate stellar masses. We than repeat our calculations with fixed photometric redshifts excluding Spitzer photometry 
and directly compare the outcomes to look for systematic deviations. Finally we use our results to compute stellar mass functions and 
mass densities up to redshift z = 5. 

Results. We find that stellar masses tend to be overestimated on average if further constraining Spitzer data are not included into the 
analysis. Whilst this trend is small up to intermediate redshifts z S 2.5 and falls within the typical error in mass, the deviation increases 
strongly for higher redshifts and reaches a maximum of a factor of three at redshift z ~ 3.5. Thus, up to intermediate redshifts, results 
for stellar mass density are in good agreement with values taken from literature calculated without additional Spitzer photometry. At 
higher redshifts, however, we find a systematic trend towards lower mass densities if Spitzer/IRAC data are included. 

Key words. Galaxies: high-redshift - Galaxies: evolution - Galaxies: fundamental parameters - Galaxies: luminosity function, mass 
function - Infrared: galaxies 



, 1. Introduction 

Whilst the assembly of stellar mass over cosmic history be- 
came a matter of interest decades ago, substantial research on 
this problem has become feasible only in the past few years. 
Extensive surveys had to be carried out, and techniques to in- 
fer stellar masses of galaxies from multicolor photometry had to 
be developed, as large complete spectroscopic samples at cos- 
mic distances are not yet available. In general, these methods 
rest upon fitting a grid of stellar population models to the data 
to compute stellar masses by multiplying the mass-to-light ra- 
tio of the best matching template with its absolute luminosity. 
This procedure has been used in numerous publications to cal- 
culate stellar mass functions or stellar mass densit i es up to in- 
termediate redshifts (e .g. Brinchmann & Enisll2000t iDrorv et alj 
120011: iBellet al.ll2003t IBorch etal4i2006) . With the increasing 
availability of deeper surveys, the analysis has been exte nded 
to higher redshifts, e.g. up to z = 2 or 3 (Dickinson et al. 2003bt 
iFontana et al.l2003l:lRudnick et al.l2003;.Fontana et al..2004). or 



even to z = 5 dProrv et alj|2005h . However, the derived results 
are typically based on observations from the U to the K-band 
and may be affected by extrapolation errors as the rest-frame 
wavelength coverage of the observed objects is shifted to the 
blue, and a good observational constraint to a galaxy's optical to 
near-infrared rest-frame luminosity is necessary to estimate its 
stellar mass reliably. Results for stellar masses may therefore be 



impaired by systematic uncertainties in the high redshift regime 
if observations are not extended to longer wavelengths. 

Since the Spitzer space telescope has become operational, it 
has been possible to complement observations up to the K filter 
with high -quality infrared data in the adjacent wavelength range 
dWerner et al. 2004). To benefit from this improvement we used 
the Great Observatories Origins Deep Survey, which provides 
deep publicly available obser vations from nurnerous f acilities in 
different wavelength regimes jDickinson et alj|2003al) . Based on 
these data we probe the influence of Spitzer photometry on the 
estimated stellar masses and examine the consequences of the 
results on inferences about the mass assembly history of the uni- 
verse. In this work we address the effects of Spitzer data to the 
mass calculation process only, i.e. we do not study whether ad- 
ditional systematic deviations of photometric redshifts emerge 
that, in general, can influence the result as well. Furthermore we 
adopt a specific set of template models for the calculations. 

The paper is organized as follows. In Sect. 2 we give a short 
overview of the dataset, and discuss the procedure adopted to 
estimate stellar masses in Sect. 3. We focus on the influence of 
Spitzer data on the derived masses in Sect. 4 and use the out- 
come to calculate stellar mass functions (Sect. 5) and mass den- 
sities (Sect. 6). Finally, we summarise our results and draw our 
conclusions in Sect. 7. 
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Throughout the paper we assume Qm = 0.3, Qa = 0.7 and 
Ho = 70 km s"' Mpc"' . Magnitudes are given in the AB system. 
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2. The dataset 

To study the assembly of stellar mass we focus on the Chandra 
Deep Field South, where numerous observations within the 
framework of the Great Observatories Origins Deep Survey 
(GOODS) provide data over a wide range of wavelengths. The 
present paper is based on GOODS- MUSIC, a multicolor catalog 
published bv iGrazian et alj (l2006l) . We briefly discuss its main 
characteristics here and refer the reader to that paper for a more 
detailed description. 

The catalog combines U-band data from the 2.2m- 
MPE/ESO telescope and VLT/VIMOS (U35, U38 and Uvim), 
Hubble/ACS images in F435W (B), F606W (V), F775W (i) and 
F850LP (z), VLT/ISAAC data in J, H and Kg, and Spitzer/IRAC 
data at 3.6 fim, 4.5 yum, 5.8 /im and 8.0 fim. Since observations 
have not yet been finished in all bands, the coverage fraction lies 
at around 63 % for UviM-data and at about 54 % for the H- 
band (ISAAC data release 1.0). Source detection has been per- 
formed independently in both the deep z-image (14651 objects 
detected) and the shallower Ks-band data (293 1 sources), end- 
ing up with a z- and Ks-complete catalog consisting of 14 847 
objec ts in total. S pecial software was developed for photom- 
etry (iDe Santis et al. 2006) to meet the requirements of color 
measurement in combined ground and space based observations 
with dissimilar point spread functions. Out of all objects, 13767, 
12041, 6767 and 5869 sources could be detected in the IRAC 
data in channels 1, 2, 3 and 4, respectively. For objects unde- 
tected in a specific image, upper limits in flux were calculated 
on the basis of morphological information derived from the de- 
tection image. As a last step redshift information were added 
to complete the dataset. To do this, spectroscopic surveys avail- 
able at that time were used to assign 1068 spectroscopic red- 
shifts. For the remaining objects, a standard photometric red- 
shift code was applied that was able to reproduce the spectro- 
scopic redshifts with an acciiracy of (1^1) - 0.045. As shown 
in Fig. 12 of iGrazian et al.l (l2006h . less than 2 % of the spec- 
troscopically observed objects reveal a photometric redshift that 
deviates severely {\zspec - Zphotl > 0.3 - 5 ■ (cr)) from the spec- 
troscopic one. 

In summary, the catalog used here consists of 14 847 ob- 
jects enclosing at least 72 stars and 68 AGNs over a total area 
of 143.2 arcmin^, with mean limiting magnitudes of ^ 26.0 
and Ks um ~ 23.8 at 90 % completeness level. 

3. Deriving stellar masses 

To calcula te galaxy st e llar ma sses we adopted the method de- 
scribed in iDrorv e t al. (2004a) and locally tested against spec- 
troscopic results in Drory et al. (2004b). This method is based 
on comparing object colors to those of a template hbrary of 
stellar population synthesis models. The five-dimensional model 
grid used here was computed from synthetic Bruzual & Chariot 
(1200 3) models with an underlying Salpeter initial mass function 
(IMF) truncated at 0.1 and 100 Mq. It is parameterized by a star 
formation history (SFH) of the form iff(t) oc exp(-f/T) evalu- 
ated at T = {0.5, 1.0, 2.0, 3.0, 5.0, 8.0, 20.0) Gyr at 15 different 
ages t € [0.2, 10.0] Gyr, with dus t extinction be tween Ay - 
and Ay -1.5 magnitudes using a ICalzetti et al] (I2OOOI) extinc- 
tion law. While covering the physical relevant parameter range 
in t and r, the upper limit in Ay must be considered as a re- 
striction. We adopted this to take into account the degeneracy in 
age and extinction to suppress solutions with unexpectedly high 
values for Ay. However, our sample may contain a small num- 
ber of heavily dust enshrouded galaxies that, in turn, will not be 



treated appropriately. In addition to the main component, a star- 
burst was superimposed which was allowed to contribute at most 
20 % to the z-band luminosity in rest-frame. It was modeled as 
a 50 Myr old episode of constant star formation with an indepen- 
dent extinction up to Ay{b) = 2.0 magnitudes. To take into ac- 
count polycyclic aromatic hydrocarbon (PAH) emission, which 
becomes important at rest-frame wavelength /I ^ 6 fim and is at- 
tributed to starforming regions, the spectral energy distribution 
(SED) of the burst compone nt was modified by including PAH- 
emission features following ,Dopita et al.l (l2005b . We found that 
the object fluxes are reconstructed slightly better at low redshifts 
compared to a burst-SED lacking this feature, while the actual 
result for calculated stellar masses is not significantly aff'ected 
(|log MwithPAH/MwithoutPAHl < 0.003 dex at z < 1). Because 
of the well known age-metallicity degeneracy, we restricted our 
models to solar metallicities and performed tests to ensure that 
we do not introduce significant systematic deviations with this 
constraint. 

To derive stellar masses, we used photometry in the filters 
U38, B, V, i, z, J, H, Ks and the four IRAC channels. In the 
blue we focused on U38 only, because the U35-filter is known 
to be lea king and the Uvim o bservations do not cover the whole 
field (see IGrazian et al.ll2006i) . We want to emphasize here that 
the Spitzer observations provide an unprecedented opportunity 
to include high-quality infrared photometry longward of about 
3 fim in stellar mass estimates. For each object we computed 
the full likelihood distribution of our models shifted to the cor- 
responding redshift. To infer the most probable mass-to-light 
ratio (M/Lk) we weighted the individual M/L^ ratios of our 
templates by their likelihoods and averaged over all parameter 
combinations. Moreover, we were able to derive an estimate of 
the expected error of this quantity from the width of the dis- 
tribution. By utilizing the Ks-band M/L ratio we benefit from 
several advantages. In general, the variation with age in M/L}^ 
is small compared to its optical counterpart as shown in Fig. [T] 
Furthermore, dust absorption is small at longer wavelengths, so 
potentially large uncertainties in Ay do not alter the result sub- 
stantially. The actual mass, which was calculated by multiply- 
ing the M/Lii ratio with the k-corrected total Ks-band lumi- 
nosity of the best fitting SED model, turns out to be compara- 
tively robust with a mean error of criogiyi = 0.16 dex. Besides 
the error attributed to the template fitting process only, calcu- 
lated stellar masses are affected by further important sources of 
uncertainty. To take them into account we performed extensive 
Monte Carlo simulations. We analyzed 1000 simulated catalogs 
where we have considered errors in the object redshifts, calcu- 
lated M/LK-ratios and a general uncertainty in photometry. The 
errors were propagated to the results for stellar masses and the 
mean uncertainty was estimated to be about criogM - 0.33 dex. 

In Fig.|2]we show the distribution of all galaxies in the mass 
versus redshift plane. We found five galaxies at high redshifts 
z > 4 with masses log M/Mq > 11.5. The presence of such 
high-mass objects less than 2 Gyr after the big bang is still con- 
troversial (see, for example. iDunlop et al.|[2007 , and references 
therein). However, a closer analysis reveals that they are more 
likely to be heavily dust enshrouded galaxies at intermediate 
redshifts. While keeping the original upper limit in Ay for the 
remaining sample, to account for the degeneracy in age and ex- 
tinction as explained above, we extended the parameter space by 
allowing a larger maximal dust extinction of Ay - 4 magnitudes 
for these five galaxies. Refitting the objects' redshifts with our 
full template library leads to much smaller redshifts z ~ 2 and 
typical masses about log M/Mq w 10.5 for four of the galaxies. 
With this reanalysis, the fit quality of our best matching SED 
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Fig. 1. Mass-to-light ratio of our template models in solar units. 
The evolution of M/Lk for different SFH r = {0.5, 1.0, 2.0, 3.0, 
5.0, 8.0, 20.0) Gyr (solid lines, from top to bottom) is shown. For 
comparison, the much wider dynamical range of the correspond- 
ing M/Ly ratio is also indicated (area between dashed lines). 
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Fig. 2. Stellar masses as a function of redshift. The typical error 
in mass is indicated in the lower right-hand corner. 



template increases significantly and becomes comparable to the 
mean for these four objects. The fifth galaxy remains at high 
z primarily because of an unexpected low Ks-band luminosity, 
which may be a result of incorrect photometry. Nevertheless, we 
want to emphasize that it is not possible to draw an unambigu- 
ous conclusion about the nature of these objects on the basis 
of currently available data. This uncertainty reflects the general 
character of samples with mainly photometric redshifts. While 
they can provide accurate data within a statistical analysis, for a 
dedicated study of rare object properties additional information 
is still required. To restrict our sample to a uniform parameter 
space, we excluded the five questionable galaxies from our anal- 
ysis in Sect. 4. 



4. The role of Spitzer 

To study the impact of Spitzer/IRAC data on our analysis, we 
considered only those galaxies with errors in mass from the fit- 
ting process criogiyi < 0.2 (^ 13 000). With this restriction, about 
half of the objects were detected in at least three IRAC chan- 
nels. However, even if a source has not been identified in a spe- 
cific filter, an upper limit in flux was included into the analysis, 
which provides a valuable additional constraint and in general 
affects the derived stellar mass. For this sample we repeated the 
calculation of stellar masses as described in Sect. 3 but reduced 
photometric information to the filters U38, B, V, i, z, J, H and Ks 
only, i.e. excluded the four IRAC channels. This subset of our 
data represents a multi-band catalog with wavelength coverage 
from about 3500 A to 23 000 A, which was typically available 
to study the high redshift universe before facilities such as the 
Spitzer space telescope became operational. We were then able 
to study the influence of additional infrared photometry by com- 
paring the properties of the best-fitting models directly. 

We found stellar masses to be overestimated without Spitzer 
data on average. The mean deviation between the masses cal- 
culated with all filters up to IRAC channel 4 (Mu-4) and 
those with limited photometry (Mu k) of the whole sample is 
log Mu-4/Mu-K - -0.18 dex with a r.m.s. of 0.41. As shown 
in Fig. |3] this trend is relatively independent of mass itself, but 
has a strong dependency on redshift. Whilst the mean deviation 
up to a redshift of z 2 is only moderate (| log Mu-4/Mu-kI < 
0.2 dex) and comparable to the typical error in mass, the dif- 
ference reaches a maximum of | log Mu_4/Mu-kI ~ 0.5 dex at 
z ~ 3.5. This corresponds to an overestimation of mass of more 
than a factor of three on average, whereas the scatter for in- 
dividual sources is large in general, as indicated by the rm.s. 
of the di stribution. Therefore we can confirm the conjecture 
made by iFontana et al.l (l2006l) . who did the same analysis on 
the Ks-selected subsample with about 3000 galaxies, although 
we found the effect of Spitzer/IRAC data to be somewhat more 
pronounced. At higher redshifts the difference in mass decreases 
again. We can therefore conclude directly that without further 
data points the extrapolation of object luminosities to longer 
wavelengths is not straight forward and induces systematic devi- 
ations. This result has been derived for a specific set of template 
SEDs with solar metallicity, a Calzetti extinction law, Salpeter 
IMF, and for fixed photometric redshifts. Although other choices 
of model parameters may result in different absolute values for 
stellar masses, we expect the detection of a systematic trend 
to be comparatively robust, as it depends on a relative devia- 
tion. T his statement is supported by the results of iFontana et al.l 
(12006 ') who found a similar trend despite using a different model 
template library and SED fitting procedure. 

The immediate reasons for the deviation can be seen in 
Fig.m where shifts in both the underlying Ks-band luminosities 
and the mass-to-light ratios are evident. For a robust calculation 
of stellar masses with the aid of infrared luminosities, a good 
knowledge of the slope of the spectra beyond the 4000 A break 
is mandatory. Whereas the brightness of an old massive galaxy 
decreases slowly in this wavelength range, a young low-mass 
object has a steep spectrum and the decline is much stronger 
To distinguish between early and late type galaxies solely on the 
basis of the blue part of their spectra is diflicult, because the pos- 
sible presence of a starburst can alter a galaxy's UV luminosity 
substantially. Furthermore, uncertainties in extinction strongly 
affect this wavelength range, which makes it even more diffi- 
cult to reveal the fundamental properties of the underlying SED. 
Therefore several data points longward of 4000 A are required 
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Fig. 3. This figure shows the difference in masses calculated un- 
der the inclusion of the Spitzer bands (Mu-4) and with restricted 
photometry (Mu-k) as a function of redshift. Averaged over all 
redshifts, the dependency on stellar mass itself is also indicated 
(inset). Mean values (solid lines) and r.m.s. of the data-points 
(dashed lines) are shown. 




Fig. 4. Ratio of inferred absolute Ks-band luminosities (upper 
panel) and shift in corresponding mass-to-light ratios in solar 
units (lower panel) as a function of redshift. The mean values 
(solid lines) and the rm.s. of the data-points (dashed lines) are 
shown. 
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Fig. 5. Correlation between the shift in masses (x-axis) and that 
in ages (y-axis) after inclusion of Spitzer photometry. The fre- 
quency of a specific combination is indicated by the color in 
nonlinear scaling. The correlation coefficient of the two quanti- 
ties is r = 0.6. 



for a reliable estimation of stellar masses. However, at a redshift 
of z = 2.5 the 4000 A break Hes at 14000 A, i.e. it has already 
passed the J filter, so the calculation is resting upon H and Ks- 
band photometry predominantly, as the slope of the rest-frame 
optical/near-infrared part of a galaxy's spectrum is characterized 
by the H - Ks color only. At a redshift of z ~ 3.5, where the de- 
viation is largest, the spectrum longward of the 4000 A break 
in rest-frame is covered solely by the Ks-filter and is therefore 
insufficiently sampled. In general, photometry becomes increas- 
ingly uncertain at higher redshifts when objects fade, so the con- 
straints on the slope weaken further. At this redshift, the sys- 
tematic deviation in masses due to the inclusion of Spitzer data 
starts to dominate over possible intrinsic errors. This effect can 
be traced back to the assignment of model SEDs that are too old 
with too high an absolute infrared luminosity. A decrease in stel- 
lar mass is therefore often associated with attributing a younger 
model SED, i.e. there is a correlation between shift in masses 
and ages, as can be seen in Fig. |5] Since extinction influences 
the shape of the spectrum in a similar way to age, the degener- 
acy in the two quantities reduces the correlation. We show the 
change of the best fitting model SED for several galaxies in the 
appendix. 

The reasons for decreasing differences in stellar masses for 
even higher redshifts beyond z = 3.5 are twofold. First, the max- 
imal age of a possible fit model is restricted by the age of the 
universe at that redshift, which acts as an upper limit. Therefore, 
the fit models are forced to be younger and a large decline in age 
with a lower resulting mass is therefore not likely. Secondly, at 
these extreme redshifts sources are very faint and an increasing 
fraction of our sample is hardly detected in the shallower J, H, 
and Ks-band data. In this case, the calculation of masses with- 
out Spitzer data is mainly based on the V, i and z photometry 
(as galaxies are U and B-filter dropouts) and it turns out that 
very young low-mass model SEDs were assigned. With the in- 
clusion of Spitzer bands the masses for this type of objects tend 
to increase, and the average difference of the whole sample gets 
smaller We want to mention here that whereas the former effect 
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Fig. 6. Deviation in mass as a function of redshift after including 
Spitzer filters one by one. Starting from masses calculated with- 
out Spitzer photometry (upper left panel, same as Fig. the 
remaining deviation is shown for calculations based on photom- 
etry including Spitzer/IRAC data {up to channel 1, upper right; 
to channel 2, lower left; to channel 3, lower right). Mean values 
(solid lines) and the r.m.s. of the data-points (dashed lines) are 
plotted. 



is universal, the latter is a property of our specific catalog and 
may be less pronounced in other surveys. 

In the next step we examined the influence of single Spitzer 
filters on the resulting stellar mass estimates. We repeated our 
calculations after adding the four Spitzer/IRAC bands to the 
analysis one by one. The outcome is illustrated in Fig. |6] where 
we show deviations in mass in relation to those calculated with 
full photometric information. Only one further data point at 
3.6 yum can reduce the remaining shift in masses to less than 
I log Mu-4/Mu-i I < 0.25 dex over the full redshift-range. After 
additionally taking into account Spitzer channel 2 at 4.5 fim, the 
systematic deviation shrinks to |log Mu-4/Mu-2l < 0.1 dex and 
thus becomes comparable to possible intrinsic errors attributed 
to the method used to derive stellar masses. Similarly, it is pos- 
sible to reduce the deviation including just Spitzer channels 3 
and 4 in the calculation. Although the remaining difference is at 
most I log Mu-4/Mu-K+3-4l ^ 0.24 dex and therefore larger than 
in the case discussed above, and taking into account lower qual- 
ity data, providing upper limits in flux can often significantly 
improve the result. Despite reduced mean deviations, the change 
in mass for single objects can still be large. 

5. The stellar mass function 

Because we found systematic shifts in calculated stellar masses, 
we wanted to examine the effects on resulting quantities such 
as the stellar mass function, i.e. the number of objects per co- 
moving volume and mass interval. To do so, we subdivided our 
catalog into seven redshift bins from z = 0.25 to z = 5. Using 
this relatively coarse grid we still work with 380 objects in the 
highest z-bin and ensure that our results are based on statisti- 
cally meaningful samples. Within each interval we calculated 
stellar mass functions after correcting our data points for incom- 
pleteness due to the flux-limited obj ect selection utilising the 
y/Vniax-formalism of ISchmid3 (1 19681) . In this context it is im- 
portant to point out that we do not achieve a well defined com- 



pleteness limit in stellar mass, as even for a sharp underlying 
flux limit, the dynamic range of mass-to-light-ratios results in 
a wide distribution in mass. Therefore the largest M/L-ratio at a 
certain redshift determines the limit in mass below which incom- 
pleteness starts to play a role. To quantify this effect we followed 
two different approaches. First, motivated by a more theoretical 
point of view, we studied the evolution of a passive evolving 
galaxy with subsolar metallicity and zero extinction formed at a 
redshift of z - 10 in a single burst. By assuming an underlying 
flux limit of ziiin = 26.0 magnitude^ we computed the corre- 
sponding stellar mass as a more conservative estimate for the 
completeness limit of the catalog. We found that galaxies with 
masses down to Mn^ w 3-10^ Mq can be detected at z - 1. 
The value is increasing to Miin, w 3 ■ 10'" Mq at z = 2.5 where 
those objects would be detected in the shallower Ks-band with 
an average flux limit of Ksum - 23.8 magnitudes. This limit in 
mass is also obtained for old but moderately star-forming sys- 
tems with extinctions around Ay * 1. In the second approach, 
which is motivated by the data actually available, we studied the 
evolution of the mass-to-light ratios of the catalog and calcu- 
lated the 95 % quantile in M/L~ as a function of redshift e.g. 
the limit below which 95 % of the M/L- ratios of our sample 
are located. Assuming a sharp flux limit such as adopted above, 
we computed the corresponding stellar mass as an estimate for 
completeness. For a redshift of z = 1 we found this value to be 
Miim w 2 ■ lO** Mq increasing to Mn^ w 9 ■ lO** M© at z = 2.5. 
Although this more aggressive method results in a lower mass 
limit, it should be clear that a significant fraction of heavily dust 
enshrouded galaxies with large M/L-ratios may stay undetected 
at intermediate z, while a reliable detection of typical more mas- 
sive active galaxies with low to intermediate extinctions should 
be possible up to high redshifts. As the two methods result in 
different completeness limits we performed the calculations in 
this section for both values independently. 

To estimate the errors of the stellar mass function we used 
1000 realisations of randomly drawn catalogs for which we have 
considered errors in computed redshifts, calculated M/Lx-ratios 
and a general uncertainty in photometry. Finally, to cope with 
incompleteness at lower masses, in view of our subsequent anal- 
ysis, we fitted our datapoints from the massive end down to the 
completeness lim it with an analytical expression suggested by 
ISchechted (Il976h of the form: 

iA(M; f, M\ a) = log(\0) ■ f ■ [lO**^-*^*>]"^"' ■ exp [-10<^-^">] 

In this formula the number density i/f(M) is parameterized via a 
scale factor ^* , a typical mass lO*' Mq and a slope-parameter 
a. First we computed the values of the fit parameters through 
an ^^-analysis in the redshift bins up to z = 4 independently. 
We excluded the highest redshift interval at 4 < z < 5.01 from 
the procedure, as we found the Schechter function to be insuffi- 
ciently constrained due to the increasing mass limit. However, a 
robust estimation of the parameters is difficult as they are highly 
degenerated in the expression used. To deal with this problem we 
decided to fix the slope-parameter at its error-weighted mean 
value, as we found no clear evidence for an evolution with red- 
shift as shown in the left-hand panel of Fig. [T] We also adopted 
this procedure because an undersampling of low mass objects in 
the catalog may affect the determination of the slope at higher 
redshifts in a systematic manner. With this constant value of a 



' The magnitude limits were calculated in lGrazian et al.l ( 2006) using 
simulations. In the z-band, the limit varies little over the field as the 
exposure map is relatively homogeneous. The variation in the Ks-filter 
is larger. 
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we repeated our calculation of the two remaining parameters M* 
and (f)* in every z-bin; for both completeness limits the result is 
plotted in Fig. [7] and listed in Table [T] A relatively uniform de- 
crease in (p* with redshift is clearly evident, which reflects the 
fact of a general decline in the number of detected objects in 
place. In contrast to the evident decrease of 0*, it is hard to say 
whether there is a hint at a mass dependent evolution of the num- 
ber density, which would manifest itself as a shift in the param- 
eter M* . Although this trend would be expected by the downsiz- 
ing scenario for galaxy evolution, where massive galaxies tend to 
form earlier than their low-mass counterparts, a slight increase 
in the typical mass with redshift followed by a decrease as indi- 
cated in the plot may also be the result of large-scale structure 
within the observed field. 

The Schechter functions can be seen in Fig. |8] where we 
show the fi t to the da t a wit h respect to the local stellar mass 
function of Cole et al.l (12001 ) with the parameters M* - 11.16, 
(/>* - 0.0031 an d a - -1.18. For com parison we also plot the 
mass function of iFontana et al.l ( l2006l) . as derived from the Ks- 
selected subsample of the GOODS-MUSIC catalog in slightly 
varied redshift bins up to z < 4. Besides a tendency to a smaller 
number of high mass objects, the datasets show a general agree- 
ment. The discrepancy at the high mass end turns out to be ro- 
bust. Although results in this region rely on only a few galax- 
ies, the error in mass is small for most of them. We discuss 
possible reasons for the dev iations in de tail in the next section. 
We also display the result of jProrv et al.l ( 120051) calculated from 
file FORS Deep Field (FDF) and a subarea of the GOODS-S 
(hereafter GSD) in the same redshift intervals without additional 
Spitzer infrared data. While the GSD data can be reproduced 
well up to intermediate redshifts, the analysis of the FDF sam- 
ple tends to result in larger values for stellar mass function. At 
high redshifts, z > 3, deviations are clearly visible for both the 
FDF and GSD datasets. They can be attributed to distinct effects; 
in addition to a shift in the mass scale as a result of the influence 
of Spitzer photometry, the total number of objects in the highest 
redshift intervals derived from an integral over the mass function 
is larger. Performing a Kolmogorov-Smirnov test on the redshift 
distributions of the FDF galaxies and the GOODS-MUSIC cat- 
alog used in this work, we can clearly reject the hypothesis that 
the two samples are drawn from the identical parent distribution 
at a 1 % level. However, both effects can have the same origin 
since the calculation of photometric redshifts becomes less reli- 
able at high redshifts if the spectra are insufficiently constrained 
in the rest-frame optical, especially as spectroscopic redshifts 
used for comparison become very rare and are restricted to the 
most luminous objects. 

6. The stellar mass density 

We are now able to compute stellar mass densities, i.e. the mass 
in stars and remnants per comoving volume at a specific redshift, 
on the basis of our results from Sect. 5. For this calculation we 
divided our stellar mass functions in each redshift bin into two 
mass intervals at the threshold values of the two completeness 
limits considered here. Above the limit, where data points and 
Schechter function fall together, we summed the stellar masses 
of our objects directly. In contrast to this procedure, we inte- 
grated the Schechter function in the low-mass range down to 
zero to take account of the fact that our catalog suffers from in- 
completeness here. The completion to lower masses contributes 
about 13 % (43 %) to the final mass density at the redshift inter- 
val 3.01 < z < 4.01 using the mass limits derived from a M/L- 
ratio analysis (passive evolution scenario). In the highest red- 



shift bin we summed stellar masses directly, not correcting for 
completeness. Although we certainly underestimate the resulting 
outcome f or stellar ma ss density, a comparison to the values pub- 
lished by iDrorv et al.l (120051) is still possible as the results were 
calculated without corrections there. In order to check our results 
for robustness, we dropped our assumption of a constant slope 
parameter a and recalculated stellar mass densities, but did not 
find appreciable deviations. To assign errors to the resulting data 
points we again used Monte Carlo simulations, and additionally 
considered uncertainties in the Schechter-parameters that affect 
the contribution of the integral over the low mass range only. 
However, a more careful analysis reveals that resulting errors 
may not include all sources of uncertainty. For example, cosmic 
variance can alter our result on a 20 % level when estimating 
the expected unce rtainty in number density of observed objects 
(ISomerville et al.l 12004). although we were able to draw on a 
relatively large survey area of about 140 arcmin- for our calcu- 
lations. Another source of uncertainty is the proper treatment of 
stars in the post-AGB phase and their infl uence on stellar popu- 
lation synthesis models (see, for example. ^ Maraston et al.ll2006l: 
Bruzual 2007; van der Wei et al. 2006). Furthermore, deviations 
from the assumed IMF can affect the outcome in a systematic 
way. In addition, it is important to point out that we can only 
give lower limits to the stellar mass densities, as we are not able 
to detect heavily dust enshrouded galaxies with large extinctions 
already at intermediate redshifts. 

We list our results in Table |2] where we also compare th e 
mass densities to the local value derived by ICole et al.l ( 1200 lb . 
It turns out that at a redshift of (z) = 1 at least 42 % of today's 
stellar mass density is already in place. This fraction decreases 
to 22 % at (z) = 2 and about 6 % at (z) = 3.5. A comparison 
with values from literature derived on the basis of a photometric 
catalog without additional Spitzer/IRAC data is shown in Fig.|9] 
Whilst up to intermediate redshifts the stellar mass density is 
reproduced well (though with much smaller scatter because of 
a larger observed area), we find systematic deviations at high 
redshifts to lower densities, which one would expect from the 
properties of the stellar mass functions discussed in Sect. 5. 

A comparison of this work with the result of iFontana et al.l 

(|2006|) . who used the Ks-selected subsample of an almost 
identical catalog with about 3000 galaxies and integrated the 
Schechter function using smoothed parameters, shows system- 
atically higher values for the stellar mass density. This trend is 
strengthening from log pthis woik/PFontana et ai. = 0.1 1 dex at a red- 
shift of (z) - 0.5 to 0.31 dex at (z) = 3.5. The discrepancy 
may have its origins in slightly different model grids used to in- 
fer mass-to-light-ratios, and in particular features of the utilized 
codes themselves . In contrast to the proceedings ofFonta na et al.l 
(I2006h . we restricted our template library to solar metallicity, 
but allowed for an independent burst component when fitting 
the object luminosities with model SEDs. As the difference in 
the inferred densities becomes more pronounced with redshift, 
it stands to reason that the derived masses of younger galax- 
ies, in particular, are subject to systematic deviations. In general, 
they are affected by larger errors as the function of the mass-to- 
light-ratio steepens at low ages. Against the background of the 
discussed degeneracies in age, metallicity and extinction, dif- 
ferences in the underlying template models can affect calculated 
stellar masses more distinctly here. Similarly, an additional burst 
component can influence the derived stellar mass. If a galaxy re- 



^ Additional spectroscopic redshift information for abo ut 150 objects 
became available in the meantime and were included in IFontana et al.l 
l l200g) . 
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Table 1. Schechter parameters of stellar mass functions as derived using the mass limits from a M/L-ratio analysis {upper section) 
and a passively evolving scenario (lower section, see text). 
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Fig. 7. Evolution of Schechter parameters of the stellar mass function with redshift, showing the slope-parameter {left panel) as de- 
rived from an independent fit in each redshift interval. The remaining parameters 0* {center) and M* {right panel) were recalculated 
using a fixed value for a. They are interdependent in that a shift in M* to higher values, for example, would result in a smaller value 
of ip*. The results using mass limits from a M/L-ratio analysis {filled circles) and a passively evolving scenario {open circles, see 
text) are plotted. 



veals both a high UV luminosity due to a recent starburst and 
a large infrared luminosity i.e. substantial mass in an old popu- 
lation, a single component fit may not be able to reproduce the 
spectrum in the whole wavelength range simultaneously, as a 
young model is too faint in the infrared and an old model not 
bright enough in the b lue. As a conseq uence, the inferred stellar 
mass can be lower (IWuvts et al.l 12007 ). Therefore, splitting up 
the fit in burst and main component covers the range of mass- 
to-light-ratios in a more flexible fashion. 



7. Summary 

In this work we estimated the influence of newly available in- 
frared data longward of the K-filter on stellar mass estimates. 
To do so we used t he GOODS-MUSIC catalog published by 
iGrazian et al.l (|2006|) . which combines photometric data in 10 
filters from 0.35 to 2.3 /urn with observations from the IRAC 
instrument of the Spitzer space telescope at 3.6, 4.5, 5.8 and 
8.0 p.m. The catalog consists of 14 847 objects within an area 
of 143.2 arcmin^ detected either in the z or the Ks-band. We 
computed stellar mass es of this sample by fitting stellar popula- 
tion synthesis models CBruzual & Charloal2003h to the data and 
multiplying the k-corrected absolute Ks-band luminosity with 
the M/LK-ratio of the best fitting model SED. To probe the in- 
fluence of the IRAC data on the analysis we repeated the compu- 
tation of stellar masses without Spitzer photometry, keeping the 
photometric redshifts fixed, and compared the outcome directly. 



Table 2. Stellar mass densities, as derived using the mass lim- 
its from a M/L-ratio analysis {upper section) and a passively 
evolving scenario {lower section, see text). The redshift interval 
4.01 < z < 5.01 has not been corrected for completeness. 



redshift interval 


logp(z)/lMQMpc-'] 


logp(z) 


P(c) 
P(z=0) 


0.25 < z < 0.75 


8.57 


0.03 


68.2 % 


0.75 < z < 1.25 


8.37 


0.02 


42.2 % 


1.25 <z< 1.75 


8.22 


0.03 


30.4 % 


1.75 <z<2.25 


8.10 


0.04 


22.9 % 


2.25 < z < 3.01 


7.93 


0.04 


15.4 % 


3.01 <z<4.01 


7.59 


0.05 


7.0% 


4.01 < z < 5.01 


>6.90 


0.08 


1.4% 


0.25 < z < 0.75 


8.57 


0.03 


68.1 % 


0.75 <z< 1.25 


8.36 


0.02 


42.2 % 


1.25 <z< 1.75 


8.22 


0.03 


30.4 % 


1.75 <z<2.25 


8.08 


0.06 


21.9 % 


2.25 < z < 3.01 


7.89 


0.06 


14.1 % 


3.01 <z<4.01 


7.61 


0.25 


7.4% 


4.01 < z < 5.01 


>6.90 


0.08 


1.4 % 



We found stellar masses to be overestimated on average, if 
further constraining infrared data from Spitzer were not included 
in the calculation. Whilst this trend is almost independent of 
mass itself, a closer analysis reveals a strong dependency on red- 
shift. While up to z w 2 the systematic deviation in mass is only 
moderate (| log Mu_4/Mu-kI < 0.2 dex) and comparable to the 
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Fig. 8. Evolution of stellar mass function with redshift, showing data points computed using the y/Vmax-formalism in seven different 
redshift intervals (filled circles) and, for z < 4, the best fitting Schechter function s utilizing mass limits from a M/L-ratio analysis 
(solid line, see text). For comparison we also plot the result of 'Drorv et al."'2005' (open triangles FDF, open squares GSD) derived 
without additional Spitzer infrared data and the result of Fontana et al. (2006), calculated within slightly different redshift bins up 
to z < 4 (open circles), and the local Schechter function of lCole et alj (120011) (dashed line). 



intrinsic uncertainty of the method adopted to estimate stellar 
masses, it increases strongly for higher redshifts and reaches a 
maximum of a factor of three at z x 3.5. The reason for this 
systematic shift can be traced back to insufficient constraints on 
the slope of the spectra redward of the 4000 A break at high red- 
shifts. It turns out that, on average, models that are too old, with 
excessively high absolute infrared luminosities and M/L-ratios 
were assigned to the data if Spitzer photometry is not included. 
Thus, a shift to lower stellar masses is likely to be correlated with 
a decreasing age of the best fitting model SED. The inclusion of 
one additional data-point longward of 3 pm can akeady reduce 
the remaining error in mass significantly. 

In the next step we used our results to calculate stellar mass 
functions in different redshift intervals utilizing the V/Vmax- 
formalism of Schmidt (1968) to correct our sample partly for 
incompleteness. To assign errors we performed extensive Monte 
Carlo simulations where we considered uncertainties in the un- 
derlying M/Lx-ratios, redshifts and a general error in photom- 
etry. Afterwards, the data points were fitted via thr ee free pa- 
ramete rs using the analytical expression suggested bv lSchechteil 
( Il976l) . We found a pronounced general decrease in numberden- 
sity of all objects with redshift. Beyond that, it is hard to say 
whether there is also a mass dependent evolution. Although the 
change in computed Schechter parameter may support this posi- 
tion, the effect can be caused by large-scale structure as well. 

Finally, we computed stellar mass densities as a function of 
redshift. We summed the mass of our objects within each red- 
shift interval at the high mass end and integrated the Schechter 
function derived on the basis of our stellar mass functions to 
complete the result for lower masses. To estimate errors we 
again used Monte Carlo simulations, but we pointed out that fur- 
ther effects such as biases in stellar population synthesis mod- 
els may be dominating sources of uncertainty. By comparing 



the outcome to the local value of ICole et alj (1200 lb we found 
at least 42 % of the stellar mass density to be already in place at 
(z) = 1 . This value decreases to 23 % at (z) - 2 and about 7 % 
at (z) = 3.5. Therefore, up to intermediate redshifts our results 
are in good agreement with values taken from literature derived 
without additional Spitzer/IRAC data. However, at high redshifts 
a systematic deviation to lower densities is present as one would 
expect from the effect of Spitzer photometry on the calculation. 
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helped to improve the presentation of our results. We are grateful to Niv Drory 
for providing the program used here to calculate stellar masses an d for valuable 
discus sions at the final stage of this paper We further acknowledge lGrazian et alj 
120061) for making the GOODS-MUSIC catalog publicly available. 

Appendix A: Examples of SED-fits 

We show some examples of a change in the best fitting model 
SED due to inclusion of Spitzer photometry. The shift in cal- 
culated stellar masses becomes manifest in a change of age, 
SFH and extinction (Figs. [ATI |A. 21 ) or the overall normalization 
(Fig. IA.3I I. In general, even a slight variation of the Ks-band 
luminosity can change the derived stellar mass substantially if 
IRAC photometry is not taken into account. On the other hand, 
we occasionally found the slope of the spectra at longer wave- 
lengths insufficiently constrained because of large photometric 
errors (Fig. |A.4| ). 
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Fig.A.l. Example of a change in the best fitting model SED 
(solid line) due to inclusion of Spitzer photometry. Upper panel: 
Without Spitzer/IRAC data the observed flux (filled circles) 
is fitted with a young starburst (dashed line) superimposed 
to a 1.5-Gyr-old main component with large extinction (dot- 
dashed line). Lower panel: After factoring in complementing 
infrared data, extinction decreases and the best fitting model 
SED becomes younger The resulting stellar mass reduces from 
log Mu-k/M0 = 10.32 dex to log MU-4/M0 = 9.69 dex. 
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Fig.A.2. Same as Fig. lA.ll The change in the fit to photo- 
metric data where the extinction varies is shown. The result- 
ing stellar mass reduces from log M\j-k/Mq - 10.07 dex to 
log Mu-4/Mo = 9.54 dex. 
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Fig.A.3. Same as Fig. lA.ll The change in the fit to pho- 
tometric data where the absolute normalisation changes dis- 
tinctly is shown. The resulting stellar mass decreases from 
log Mu-k/M© = 9.26 dex to log Mu-a/Mq = 8.52 dex. 
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Fig.A.4. Same as Fig. lA.ll The figure shows an example of an 
insufficiently constrained optical- to near-infrared-luminosity 
due to large photometric errors. The resulting stellar mass de- 
creases from log Mu-k/Mq - 9.45 dex to log Mjj^^/Mq - 
8.65 dex. 



